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Dynamical chiral-symmetry breaking at T = and T^O 
in the Schwinger-Dyson equation with lattice QCD data 
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For the study of dynamical chiral-symmetry breaking (DCSB) in QCD, we investigate the 
Schwinger-Dyson (SD) formalism based on lattice QCD data. From the quenched lattice data 
for the quark propagator in the Landau gauge, we extract the SD kernel function K(p 2 ), which is 
the product of the quark-gluon vertex and the polarization factor in the gluon propagator, in an 
Ansatz-independent manner. We find that the SD kernel function K(p 2 ) exhibits infrared vanishing 
and a large enhancement at the intermediate-energy region around p ~ 0.6GeV. We investigate 
the relation between the SD kernel and the quark mass function, and find that the infrared and 
intermediate energy region (0.35GeV < p < 1.5GeV) would be relevant for DCSB. We apply the 
lattice-QCD-based SD equation to thermal QCD, and calculate the quark mass function at the finite 
temperature. Spontaneously broken chiral symmetry is found to be restored at the high temperature 
above 100 MeV. 

PACS numbers: 12.38.Aw, 12.38Lg, 12.38Mh, 12.38-t. 



I. INTRODUCTION 

Quantum chromodynamics (QCD) has been accepted 
as the fundamental theory of the strong interaction of 
hadrons Due to the asymptotic freedom of QCD, 
the perturbative calculation is applicable to the high- 
energy process in the hadron reactions. In contrast, 
low-energy QCD becomes a strong-coupling gauge the- 
ory and exhibits interesting nonpcrturbative phenom- 
ena such as color confinement 13 and dynamical chiral- 
symmetry breaking (DCSB) 

In perturbative QCD, the QCD scale parameter Aqcd 
is one of the most important quantities, and provides a 
typical scale of the strong interaction. At the one-loop 
level, the running coupling constant is expressed as 



9 2 (P 2 

4-7T 



12tt 



{UN c -2N f )Hp 2 /A 2 QCD ) 



(1) 



where N c and Nf being the numbers of the color and 
flavor, respectively 0. 

The perturbative result is valid only at a large mo- 
mentum scale p 2 . Indeed, the running coupling constant 
a s (p 2 ) becomes large in the infrared region, and formally 
diverges at p 2 = Aq CD in Eq. . Accordingly, the simple 
perturbation is no more meaningful and some nonpcrtur- 
bative effect should appear in this strong-coupling region. 

As a quantitative aspect, the experimental data sug- 
gest a s {M%) ~ 0.1172 ± 0.002 and the QCD scale pa- 
rameter is estimated as Aqcd (MS) — 216 ± 25MeV in 
the MS scheme with Nf = 5 8}. From Eq.JTJ, one finds 
a rough estimation as a s ((0.5GeV) 2 ) ~ 1, which means 
the breakdown of the validity of the perturbation the- 
ory at p ~ 0.5GeV. In other words, one may expect the 
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appearance of nonpcrturbative effects due to the strong 
coupling in the infrared region as p ~ 0.5GeV 

DCSB is one of the outstanding nonperturbative fea- 
tures in QCD. The chiral symmetry which the QCD la- 
grangian possesses in the massless quark limit, is sponta- 
neously broken in the nonperturbative QCD vacuum [jj . 
In the mathematical viewpoint, DCSB leads to the low- 
energy theorem and the current algebra and, in 
the phenomenological viewpoint, it results in the absence 
of the parity doubling for hadrons. In terms of the non- 
trivial QCD vacuum, DCSB is characterized by the quark 
condensate {qq) ~ -(225±25MeV) 3 which is caused 
by attractive interaction acting the quark-antiquark pair 
like Cooper-pair condensation in the superconductiv- 
ity [3j. Such a physical process has been demonstrated 
by using the effective models of QCD such as the Nambu- 
Jona-Lasinio model 0, ^3 an d the instanton model ^| • 
These models suggest that the almost massless quark ac- 
quires a large effective mass of M ~ 300MeV as a result of 
DCSB and behaves as a massive constituent quark in 
the infrared region. The pion is identified as the Nambu- 
Goldstone boson associated with DCSB and obeys the 
low-energy theorem and the current algebra [9lll0|. where 
the pion decay constant, f n ~ 93MeV, is also a relevant 
quantity characterizing DCSB In this way, DCSB 
is characterized by several quantities, the quark conden- 
sate {qq) ~ — (225 ± 25MeV) 3 , the effective quark mass 
M ~ 300MeV and the pion decay constant f % ~ 93MeV. 

As the study of DCSB based on QCD, the Schwinger- 
Dyson (SD) equation 0. U and the Bethe-Salpeter 
equation which are expressed as nonlinear integral 
equations, have been used as popular methods to incor- 
porate the infinite-order effects on the gauge coupling. 
In most studies, however, the free gluon propagator was 
used in the SD equation, and possible nonperturbative 
effects in the infrared region were neglected. In several 
studies, the infrared nonperturbative effect was modeled 
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and was taken into account in the SD equation |l5l Il6| , 
but such modeling was done hypothetically. Recently, 
the nonperturbative quark and gluon propagators have 
been calculated in lattice QCD Monte Carlo simulations 
with high accuracy. In this paper, we try to combine the 
SD approach and the lattice QCD result for the quark 
and gluon propagators as the first step to understand 
the hadron physics in terms of quarks and gluons. 

In section 2, we summarize the quark and the gluon 
propagators in lattice QCD. Using the lattice quark prop- 
agator, we calculate the pion decay constant f v in the 
Pagels-Stokar approximation, and the quark condensate 
(qq). In section 3, we briefly review the SD formalism 
for the quark field in QCD. In section 4, from the lat- 
tice QCD data of the quark propagator, we extract the 
SD kernel function, which is expressed by the product of 
the quark-gluon vertex and the polarization factor in the 
gluon propagator. In section 5, we investigate the rela- 
tion between the SD kernel and the quark mass function 
through several tests with modified SD kernels. In sec- 
tion 6, we apply the lattice-QCD based SD equation to 
the finite temperature system, and calculate the quark 
mass function and the quark condensate. Section 7 is 
devoted to the summary and concluding remarks. 



II. PROPAGATORS OF QUARKS AND 
GLUONS IN LATTICE QCD 

In this section, we summarize the recent lattice QCD 
results for the quark and the gluon propagators. The 
lattice QCD Monte Carlo simulation is the first princi- 
ple calculation of the strong interaction directly based 
on QCD in the Euclidean metric [13. In these years, 
relatively accurate lattice QCD calculations have been 
performed for the quark and gluon propagators in the 
Landau gauge at the quenched level |lSt liH l2*cj . These 
propagators are considered to include all the nonpertur- 
bative effects in quenched QCD. 

In the Landau gauge, the Euclidean gluon propagator 
is generally expressed by 



dip 2 ) 



P 



P^Pv 
p 2 



(2) 



where we refer to d(p 2 ) as the "polarization factor" in 
the gluon propagator. 

We find that the lattice QCD data 0] for the po- 
larization factor d{p 2 ) is well described by the analytic 
function of 



d( P 2 ) = Z 



4 1 2 
p + ap 



b 2 



(3) 



with a ~ 1.94GeV 2 and b ~ 0.486GeV 2 , as shown in 
Fig.l. The constant factor Z ~ 1.80 denotes the wave- 
function renormalization of the gluon field at the lattice 
cutoff. 

Note here that there are two remarkable features in the 
functional shape of d(p 2 ). 



1. The polarization factor d(p 2 ) exhibits the infrared 
vanishing and dip 2 ) is proportional to p 2 in the 
infrared region as p < 0.5GeV. 

2. The polarization factor d(p 2 ) exhibits a large en- 
hancement in the intermediate-energy region as 
p ~ lGeV. 

Next, we summarize the quark propagator S(p) in the 
Landau gauge in lattice QCD. In general, the quark prop- 
agator in the Landau gauge is expressed as 



S(p) 



Zip 2 ) 



p + Mip 2 ) 



(4) 



in the Euclidean metric. Here, M(p 2 ) is called as the 
quark mass function, and Zip 2 ) corresponds to the wave- 
function renormalization of the quark field 0, 01 • In 
terms of the quark propagator, DCSB is characterized 
by the mass generation as M(p 2 ) 7^ 0. 

The quark mass function Mip 2 ) in the Landau gauge 
is recently measured in lattice QCD at the quenched 
level 0| . Figure 2 shows the quark mass function Mip 2 ) 
in the chiral limit obtained from lattice QCD at the 
quenched level [lq. The lattice data in the chiral limit 
can be fitted by 



Mip 2 



1 + ip/pp 



(5) 



with M ~260 MeV, p ~870MeV and 7 -3.04 [lj, al- 
though there is some ambiguity on the fit function. The 
infrared quark mass M(0) = Mq ~ 260MeV seems con- 
sistent with the constituent quark mass in the quark 
model [T^ |. 

As for the quark wave-function renormalization factor 
Zip 2 ), the lattice QCD results seem somehow unclear 0, 
l20l | due to the relatively large statistical and systematic 
errors. For large p 2 , one finds Zip 2 ) ~ 1 in the Landau 
gauge, while lattice data seem to indicate Zip 2 ) < 1 for 
small p 2 . For the practical calculation in this paper, we 
adopt Zip 2 ) ~ 1 in the Landau gauge. 

Using the lattice QCD result for the quark propagator, 
the quark condensate (qq) , which is the order parameter 
of DCSB, can be estimated from the trace of the quark 
propagator as 



(qq)a = - 



P 



Air 2 



(2tt)< 

A 2 



trSgip) 



7^ dfV 



2 2 Zip 2 )Mip 2 ) 



P 



M 2 ip 2 



(0) 



at the renormalization point A. With the assump- 
tion of Zip 2 ) = 1, one obtains (?<z)a=igcV — 
-(223MeV) 3 , (qq) A =4GcV - -(300MeV) 3 and 
(qq) a=oo — — (320MeV) 3 , which seem consistent 
with the standard value (qq) ~ -(225 ± 25MeV) 3 . 
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We calculate the pion decay constant f^, which also 
characterizes DCSB, in the Pagels-Stokar approxima- 
tion H3, 



4tt 2 



dp 2 



P 2 M(p 2 ) 



{p 2 + M 2 {p 2 )} 2 

2 dM(p 2 ) 



in the formalism jl^].) Then, one obtains 

M(p 2 ) 2 f A Z(q 2 )M(q 2 ) 

Z(p 2 ) 9 J (2Tr)U 2 + M 2 {q 2 ) 

X T((p-q) 2 )D^((p-q) 2 ). 



(11) 



M(p 2 ) 



El 

2 



dp 2 



(7) 



In the Landau gauge, the Euclidean gluon propagator 
takes the general form of 



where Z(p 2 ) = 1 is assumed. With the lattice data for 
M(p 2 ) in Eq.JSJ), we obtain /„. ~ 87MeV, which is con- 
sistent with the experimental value /„• ~ 93MeV or the 
theoretical estimation of /„• ~ 87MeV in the chiral limit 
from the chiral perturbation [TT| . 



III. THE SCHWINGER-DYSON FORMALISM 
FOR QUARKS 

The Schwinger-Dyson (SD) equation for the quark 
propagator S(p) is described with the nonperturbative 
gluon propagator D^fj)) and the nonperturbative quark- 
gluon vertex gT v {p, q) as 

S-\p) = Sv 1 + C F g 2 J A«^S(q)D^(p - q)T v {p, q) 

in the Euclidean metric. Here, So(p) denotes the bare 
quark propagator, and the color factor of quarks has been 
calculated as Cf =4/3. 

In the practical calculation for QCD, however, the 
SD formalism is drastically truncated: the perturbative 
gluon propagator and the one-loop running coupling are 
used instead of the nonperturbative quantities in the orig- 
inal formalism. This simplification seems rather danger- 
ous because some of nonperturbative-QCD effects are ne- 
glected. 

We formulate the SD equation for quarks in the chiral 
limit in the Landau gauge. By taking the trace of Eq. (JHJ), 
one finds 



d,Ap 2 ) = 'X 1 ( *„, 



(12) 



with the gluon polarization factor d(p 2 ). Therefore, 
Eq.JTTJ is expressed as 



Z(p 2 



d A q Z{q 2 )M(q 2 ) 
(2tt) 4 q 2 + M 2 (q 2 ) 
r((p-q) 2 )d((p-q) 2 ) 



(p - q) 2 

Here, we define the kernel function 

K{p 2 ) ee g 2 T(p 2 )d(p 2 ) 



(13) 



(14) 



(8) 



as the product of the quark-gluon vertex T(p 2 ) and the 
gluon polarization factor d(p 2 ). Then, the SD equation 
is rewritten as 



M{p 2 



3C f 



d*q Z(q 2 )M(q 2 )K((p-q) 2 ) 



(15) 



Z(p 2 ) r J (2tt) 4 q 2 + M 2 (q 2 ) (p - q) 2 

Note that the actual kernel Ksb(p 2 ) of the SD equation 
includes the Coulomb-propagator factor 1/p 2 as 



, _ K (p 2 ) _ g 2 T(p 2 )d(p 2 ) 



P 



P 



(16) 



In the Landau gauge, from several theoretical argu- 
ments |4j, |7J, |22j , we expect that the quark wave- function 
renormalization is not so significant, and therefore we ap- 
proximate as Z(p 2 ) = 1. In this approximation, the SD 
equation reduces to 



M(p 2 ) _ C F g 2 f d 4 q 



M(p 2 ) = 3C P 



x tr^ 7Al 



4 J (2tt) 4 p 
Z(q 2 " 



d±q M(q 2 ) K((p-q) 2 ) 
(2tt) 4 q 2 + M 2 {q 2 ) (p - q) 2 ' 



(17) 



A + M(q 2 ) 



For the quark-gluon vertex, we assume the chiral- 
preserving vector-type vertex, 



r^(p,«) = 7 M r((p- g) 2 ), 



(10) 



which keeps the chiral symmetry properly. Note here 
that, to preserve the Ward-Takahashi identity for the ax- 
ial vector vertex, the gluon momentum (p — q) should be 
taken as the argument of the quark-gluon vertex T in the 
ladder approximation of the SD and BS equations |22| . 
(In contrast, to be strict, the Higashijima-Miransky ap- 
proximation 0, explicitly breaks the chiral symmetry 



In the following sections, we analyze this type of the SD 
equation. 



IV. EXTRACTION OF THE KERNEL 
FUNCTION IN THE SD EQUATION FROM 
LATTICE QCD 

In this section, we extract the kernel function K(p 2 ) = 
g 2 T(p 2 )d(p 2 ) in the SD equation l|17|) from the quark 
mass function M(p 2 ) obtained in lattice QCD in an 
Ansatz-independent manner. By shifting the integral 
variable from q to q = p — q, we rewrite Ea. (|17l) as 



M(p 2 



3C f 



d 4 q 



M((p-q) 2 ) 



K{q 2 \ 



{2ttY (p-q) 2 +M 2 {{p-q) 2 ) q 2 



(18) 
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Therefore, we obtain 



where 0(p , q ) is defined with M(p ) as 



M(p 2 ) = / dq 2 e(p 2 ,q 2 )K(q 2 ), (19) 
Jo 



e(pV) 



3C F 
8tt 3 



d6» sin 2 6>- 



M (p 2 + g 2 - 2pg cos ( 



p* + q 2 — 2pq cos 9 + M 2 (p 2 + q 2 — 2pq cos t 
I 



(20) 



Regarding the momentum squared (p 2 , q 2 ) as suffices (m, 
n), Ea. ljT9l can be rewritten as 

M ro = ^e TO „isr„. (21) 

n 

Here, 6 m „ is a real symmetric matrix on m and n as 

8mn = B nm S R, (22) 

because of Q(p 2 ,q 2 ) = 0(q 2 ,p 2 ) £ R. 

Once the quark mass function M(p 2 ) is obtained, 
Q(p 2 ,q 2 ) is calculable with Eq.(J20|), and, using Eci. (|2~TJl . 

we can extract K n directly from mn and M n as 



(23) 



Since the quark mass function M(p 2 ) is given by Eq.(|3J) in 
lattice QCD, we can calculate the kernel function K(p 2 ) 
from Ea. (|23[l without any assumption of the functional 
form on the kernel function K(p 2 ). For the practical 
calculation, we discretize the momentum squared p 2 and 
q 2 after the proper transformation as p 2 — arctana and 
q 2 = arctan/3, and solve Ea. (|2*T|) for K n . 

As shown in Fig. 3, we numerically extract the kernel 
function K(p 2 ) = g 2 T(p 2 )d(p 2 ) from the lattice QCD 
result, Eq.(|SJ|, for the quark propagator in the Landau 
gauge. For the check of the validity, we have confirmed 
that the mass function M{p 2 ) in Eq.(0 is precisely re- 
produced with the obtained kernel function K(p 2 ). 

As remarkable features, we find "infrared vanishing" 
and "intermediate enhancement" in the kernel function 
K(p 2 ) in the SD equation pf : 

1. The SD kernel function K(p 2 ) seems consistent 
with zero in the very infrared region as 



K (p 2 < O.lGeV 2 ) ~ 0. 



(24) 



2. The SD kernel function K(p 2 ) exhibits a large 
enhancement in the intermediate-energy region 
around p ~ 0.6GeV. In fact, K(p 2 ) takes the max- 



imal Value -ftfmax 

(0.607GeV) 2 . 



38.76 at p 2 ~ 0.368G6V 2 



These tendencies of infrared vanishing and interme- 
diate enhancement in the kernel function K(p 2 ) = 



g 2 T(p 2 )d(p 2 ) are qualitatively observed also in the di- 
rect lattice-QCD measurement for the polarization factor 
d(p 2 ) in the gluon propagator in the Landau gauge 0] 
as shown in Fig.l. 

Note here that the usage of the perturbative gluon 
propagator and the one-loop running coupling g(p 2 ) cor- 
responds to the perturbative SD kernel function 



K pcrt (p 2 )=g 2 (p 2 ) = 



48tt 2 



(11N C 



2N f )Hp 2 /A 2 QCB ) 



(25) 



For comparison, we add by the dotted curve in Fig. 3 the 
perturbative SD kernel function K pelt (p 2 ) at the one- 
loop level. Both in the infrared and in the intermediate 
energy regions, K pclt (p 2 ) largely differs from the present 
result K(p 2 ) based on lattice QCD. In fact, the simple 
version of the SD equation using the perturbative gluon 
propagator and the one-loop running coupling would be 
too crude for the quantitative study of QCD. 



V. RELATION BETWEEN DCSB AND THE SD 
KERNEL 

In this section, we investigate the relation between the 
quark mass function M(p 2 ) and the SD kernel function 
K{p 2 ) through several simple tests, considering the rele- 
vant momentum region for DCSB. Since the SD equation 
is a nonlinear equation, it is difficult to find out such rel- 
evant region for DCSB in a rigorous manner. In order to 
determine the relevant momentum region for DCSB, we 
calculate the quark mass functions for modified SD ker- 
nels and compare the results with the original one. We 
consider the following four modifications. 

1. An ultraviolet (UV) cutoff for the SD kernel 

2. An infrared (IR) cutoff for the SD kernel 

3. An intermediate (IM) suppression for the SD kernel 

4. A simple scaling of the SD kernel 

A. Ultraviolet cut case 

To check the relevant momentum region for DCSB, we 
investigate the ultraviolet (UV) cut case, infrared (IR) 
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cut case and intermediate (IM) suppression case for the 
SD kernel, To begin with, we investigate the UV-cut case 
for the SD equation 1)17(1 by using the UV-cut SD kernel 
function, 

K VV ( P 2 ; A uv ) = *Tf/)0(A uv - p 2 ), (26) 

instead of K(p 2 ). 

For each Auv, we solve the UV-cut SD equation and 
obtain the corresponding solution Muv(p 2 ; Auv) for the 
quark mass function. We show in Fig. 4(a) the in- 
frared quark mass Muv(0;Auv) plotted against Auv- 
For Auv > 2GeV, almost no effect is observed as 
Mxv(p 2 ;Auv) — M\jv(p 2 ',°o) — M(p 2 ), which clearly 
indicates that UV region is not important for DCSB. In 
contrast, a significant reduction of the quark mass func- 
tion Muv(p 2 ;Auv) is observed for Auv < 1.5GeV. In 
particular, for Auv < Ayy — 0.9GeV, no DCSB is ob- 
served as Muv(j? 2 ;Auv) = 0, which may suggest that 
the momentum scale of p ~ lGcV plays an important 
role for DCSB. Figure 4(b) shows the UV cut SD kernel 
function for Ayy, which is the critical value on DCSB. 

This result seems natural because the strong coupling 
nature at the infrared and intermediate energy regions 
would be essential for DCSB in QCD [HHI. 



B. Infrared cut case 

To investigate the role of the infrared (IR) region for 
DCSB, we consider also the IR cut case by using the 
IR-cut SD kernel function, 

K m (p 2 ;A m ) ee K (p 2 )6(p 2 ~ Aj R ), (27) 

instead of K(p 2 ) in Eq. l|17[l . For each Air,, we solve the 
IR-cut SD equation and obtain the corresponding so- 
lution Min(p 2 ; Air). We show in Fig. 5(a) the infrared 
quark mass Mir(0; Air) plotted against the IR cutoff 
Air, 

For Air <0.35GeV, no significant effect is observed for 
DCSB, according to the infrared vanishing of the SD ker- 
nel. In contrast, for Air, > Agj* ~ 0.53GeV, no DCSB 
is observed as Mir(p 2 ;Air) = 0. These results seem 
to indicate the relevant role of the infrared region as 
0.35GeV < p < 0.53GeV for DCSB. Figure 5(b) shows 
the IR-cut SD kernel function for AJr*, which is the crit- 
ical value on DCSB. 



C. Intermediate suppressed case 

The SD kernel function K{p 2 ) obtained from the lat- 
tice QCD data of the quark propagator indicates the 
intermediate enhancement, and takes a maximal value 
K max = 38.76 at p 2 ~ 0.368GeV 2 ~ (0.607GeV) 2 in the 
present analysis. (See Fig. 3) To investigate the role of 
the intermediate (IM) enhancement of the SD kernel for 



DCSB, we examine the IM-suppressed SD kernel func- 
tion, 

K 1M {p 2 -c) ee Min(AV), cK max ), (28) 

with a real and positive constant c. 

For each c, we calculate the quark mass function 
A/im(p 2 ; c) using the SD equation l|17|l with the modified 
kernel Kjyi{p ; c). We show in Fig. 6 (a) the infrared quark 
mass Mim(0;c) plotted against c. Of course, for c > 1, 
one finds Ki M {p 2 ]c) = K(p 2 ) and M IM (p 2 ;c) = M(p 2 ). 
As c decreases from 1, Miuip 2 ] c) rapidly decreases, and 
no DCSB is found as Mim(p 2 ', c) — for c < c crit ~ 0.58. 
Figure 6(b) shows the IM-suppressed SD kernel func- 
tion for Cent j which is the critical value on DCSB. This 
would indicate that the intermediate enhancement of the 
SD kernel in the region of 0.2GeV 2 < p 2 < 0.8GeV 2 , 
i.e., 0.4GeV < p < 0.9GeV, plays an important role for 
DCSB. 

From the above analyses in Subsect. A, B and C, the 
relevant momentum region for DCSB is considered to be 
the IR and the IM regions as 

0.35GeV < p < 1.5GeV. (29) 

D. Sensitivity of DCSB on the magnitude of the 
SD kernel 

In general, DCSB is realized when the SD kernel is 
large enough. In fact, the strong-coupling nature is ex- 
pected to be the origin of DCSB, as was first pointed out 
in the Nambu-Jona-Lasino (NJL) model where DCSB 
is realized for the strong coupling case, e.g., ffN.u.A 2 > 
ttN c for the quark version of the NJL model [j, yj [TJ • 

Finally in this section, we consider a simple scaling of 
the SD kernel with a real and positive constant A. In 
fact, we solve the modified SD equation with the scaled 
SD kernel function 

K(p 2 ;\) = XK(p 2 ), (30) 

and investigate the corresponding quark mass function 
M (p 2 ; A) . This modification can be physically inter- 
preted as the check on the sensitivity of DCSB on the 
gauge coupling magnitude, since the K(p 2 ) is propor- 
tional to g 2 . In fact, the SD-kernel scaling as Ea. (|30|l is 
equivalent to the scaling of the QCD gauge coupling g 2 
as 

gl = a 9 2 (3i) 

in the SD equation. (From Eq. (|15fl , this simple SD-kernel 
scaling can be also interpreted as the introduction of the 
constant quark wave-function renormalization Z as A = 
Z 2 , although we have ignored its effect as Z = 1.) 

For each A, we solve the SD equation 1|17|) with the 
scaled kernel function K(p 2 ;X), and obtain the corre- 
sponding quark mass function AI(p 2 ;X). We show in 
Fig. 7 the infrared quark mass M (0; A) against the scale 
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parameter A. For A > 0.8, the quark mass function 
M (p 2 ; A) is an increasing function on A for fixed value of 
p 2 . For A < A crit ^ 0.8, no DCSB occurs as M{p 2 ; A) = 0. 

This result indicates that the physical case of A = 1 
may locate near the critical case A cr it — 0.8 for DCSB. 
In accordance with the closeness to the critical case, the 
solution for the quark mass function seems rather sensi- 
tive on the overall scaling of the SD kernel near A = 1. 
This may be rather serious for the SD approach because, 
for more realistic calculation, we need more sophisticated 
treatment on the wave function renormalization, the ver- 
tex and so on. Indeed, a similar approach in Ref.p^] 
results in a smaller value for the pion decay constant f n 
from the BS equation with Z(p 2 ) taken into account. 

In this way, DCSB such as the quark mass function 
may be sensitive to the overall scaling of the SD kernel. 
Our result indicates that DCSB is near the critical point 
as A cr it — 0.8, i.e., g 2 ~ 1.25g 2 lit for the coupling constant 
from Eq. (j2U . 



VI. CHIRAL SYMMETRY AT FINITE 
TEMPERATURE 

Finally, we demonstrate a simple application of the 
lattice-QCD-based SD equation to chiral symmetry 



Here, Mt(p 2 ,w 2 ) depends only on p 2 and ui 2 , because 
of the three-dimensional rotational invariance and the 
imaginary-time reversal invariance of the system. 

Using the kernel function K(p 2 ) obtained in section IV, 
we solve the thermal SD equation lj3~5l) for Mt{p 2 ,uj 2 ). 
Figure 8 shows the numerical result for the thermal in- 
frared quark mass Mp(p 2 = 0,Wq) plotted against the 
temperature T. Chiral symmetry restoration is found 
at a critical temperature T c ~ lOOMeV, which seems 
rather small in comparison with the critical temperature 
T c = 260 - 280MeV of the quenched QCD phase transi- 
tion m. 

In this calculation, the chiral phase transition is found 
to be of the second order, and we find near T c 

M t (0,lo 2 )~(1-T/T c ) 0A7 (36) 

as the critical behavior, although more careful treatments 
would be necessary for the argument of the critical phe- 
nomena. 

Figure 9 shows the thermal quark mass function 
A/t(p 2 ,w 2 ) at various Matsubara frequencies at a low 



restoration in finite-temperature QCD, using the Mat- 
subara imaginary-time formalism. 

In this formalism, the quark field q(x, t) and the 
gluon field (x, r) at a finite temperature T obey the 
anti-periodic and periodic boundary condition in the 
imaginary-time direction, respectively, as 

g(x,r+ 1/T) = -q(x,r), 
A p (x,r+1/T) = A„(x,r), (32) 

Accordingly, the temporal momentum variable po is dis- 
cretized to be the Matsubara frequency, 

Po —>■ (2n + 1)ttT = uj n for quarks, 

Po —>■ 2nirT for gluons, (33) 

and the corresponding integration over pa becomes the 
summation over the Matsubara frequencies as 

J^E (34, 

~°° n=-oo 

in the SD equation. 

As a result, the SD equation for the thermal quark 
mass function A/t(p 2 ,w 2 ) of the Matsubara frequency 
u) n = (2n + l)irT is given as 



(35) 

I 

temperature T — 0.6T C = 60MeV and at a high tem- 
perature T = 0.9T C = 90MeV. One finds the following 
tendencies for the thermal quark mass Mr(p 2 ,o; 2 ). 

1. At fixed T, Mt(p 2 ,w 2 ) is a decreasing func- 
tion of p 2 for each Matsubara frequency n, and 
Mr(p 2 , w 2 ) decreases with w 2 for each value of p 2 . 

2. For fixed n and p 2 , Mr(p 2 ,a> 2 ) decreases with the 
temperature T. 

As an interesting dependence of u> 2 and p 2 , we find at 
each T the "covariant-like relation" [23 for the thermal 
quark mass function as 

M T (p 2 , w 2 ) ~ M T (p 2 + c4) = M T {p 2 ). (37) 

Here, p 2 = p 2 + u 2 = p 2 + {(2n+ 1)ttT} 2 corresponds to 
the four-dimensional momentum squared, p 2 = p 2 + Pq, 
and actually reduces into p 2 at T = 0. For the demon- 
stration of this relation, we define Mt(p 2 ) = Mx(p 2 , uJq) 
and compare Mxip 2 ) with A/t(p 2 ,w 2 ) for each n in 
Fig. 10. We find an approximate coincidence between 
them at each n even at temperature close to T c . 



00 

M T (p 2 ,ujl) = 3C F T J 



m— — 00 



d 3 q M T (q 2 ,i4j K((u; n -u; m ) 2 + (p-q) 2 ) 

(2^3a; 2 a + q 2 + M 2 (q 2 , Wl 2 n ) ( W „ - uj m ) 2 + (p - q) 2 ' 
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The quark condensate at finite temperature is calcu- 
lated with M t (p 2 ,ujI) as [H| 



{ q qi T)) A = -Tj2f-^trS q (p) 



(38) 



2N r , 



N 



n=-JV ,/U 

" AT->LO„ 



~2 M T (p 2 -ul,u>l) 



where AT = — |] is the UV cutoff on the Matsubara 
frequency. Note that Y satisfies (2Y + 1)ttT ~ A, and 
corresponds to the momentum cutoff A. In the limit of 
T —>■ 0, (qq(T))\ becomes the qua rk condensate (gg)A in 
Eq. © at the zero temperature J25j . We show the thermal 
quark condensate (qq(T))A plotted against the tempera- 
ture T in Fig. 11, which exhibits the critical behavior near 
T c as 



(qq(T)) A=4GcV ~ (1 - T/T c ) 



0.74 



(39) 



In this calculation, we have included only the (anti- 
periodicity in the imaginary-time direction, and have 
ignored the nontrivial thermal effects on the quark and 
gluon propagators and vertex functions. Nevertheless, we 
observe chiral symmetry restoration at high temperature 
and obtain a rough estimate of the critical temperature. 
For more consistent calculation, it would be interesting 
to use the lattice QCD results on the quark and gluon 
propagators at finite temperature. 



VII. SUMMARY AND CONCLUDING 
REMARKS 



We have investigated the relation between the quark 
mass function and the SD kernel, considering the impor- 
tant scale region for DCSB. Several checks have indicated 
that the infrared and the intermediate energy regions as 
0.35GeV < p < 1.5GeV would be relevant for DCSB. 
The "intermediate kernel enhancement" is found to be 
also important for DCSB. 

We have applied the lattice-QCD-based SD equation 
to thermal QCD, and have calculated the quark mass 
function at the finite temperature. We have found that 
spontaneously broken chiral symmetry is restored at the 
high temperature above 100 MeV. 

In this paper, we have investigated the SD equation 
in the Landau gauge at the quenched level. It is inter- 
esting to study the SD equation in other choices of the 
gauge such as the maximally Abelian gauge [2]]], which 
may reveal possible relations between DCSB and color 
confinement [l^. It is also interesting to investigate the 
SD equation based on the quark and the gluon propaga- 
tors in full QCD. As a next possible step, the application 
to finite density QCD would be an interesting subject, 
which is rather hard to be performed within the present 
lattice-QCD ability. 

In near future, we expect that the combination of lat- 
tice QCD and the nonperturbative formalism such as the 
SD and the BS equations provides us a powerful tool 
to clarify nonperturbative aspects of the hadron physics 
based on QCD. In this framework, if one obtains accurate 
lattice QCD data, one would derive more reliable results. 
Therefore, in order to establish this new framework, it is 
much desired to obtain more accurate data on the quark 
and gluon propagators as well as vertex functions in lat- 
tice QCD. 



We have investigated the Schwinger-Dyson (SD) for- 
malism based on lattice QCD data, and have studied dy- 
namical chiral-symmetry breaking (DCSB) in QCD, We 
have extracted the SD kernel function K(p 2 ), which is the 
product of the quark-gluon vertex and the polarization 
factor in the gluon propagator, in an Ansatz-independent 
manner from the quenched lattice data for the quark 
propagator in the Landau gauge. We have found that the 
SD kernel K(p 2 ) exhibits infrared vanishing and a large 
enhancement at the intermediate-energy region around 
p ~ 0.6GeV. 
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FIG. 1: The lattice QCD result taken from Ref.pt for the 
polarization factor d(q 2 ) = q 2 D(q 2 ) in the gluon propagator 
in the Landau gauge. 




FIG. 2: The lattice QCD result taken from Ref.Qjj] for the 
quark mass function M(q 2 ) in the chiral limit in the Landau 
gauge. 



u 1 1 1 1 1 1 1 1 — 1 
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FIG. 3: The kernel function in the SD equation, K(p 2 ) = 
g 2 r(p 2 )d(p 2 ), extracted from the lattice QCD result of the 
quark propagator in the Landau gauge. The calculated data 
are denoted by the square symbols, and the solid curve de- 
notes a fit function for them. The dotted curve denotes the 
perturbative SD kernel function K pert (p 2 ) for comparison. As 
remarkable features, K(p 2 ) exhibits infrared vanishing and in- 
termediate enhancement. 



11 




1 1 ■ J 1 1 

0.2 0.4 0.6 0.8 1 1.2 

p 2 [GeV 2 ] 



FIG. 4: (a) The infrared quark mass Muv(0; Auv)/Mo in 
the UV-cut SD equation plotted against the artificial UV- 
cutoff parameter Auv- Mo ~ 260MeV denotes the infrared 
quark mass. For Auv > 2GeV, almost no effect is ob- 
served. For Auv < Auv ^ 0.9GeV, no DCSB is observed as 
M uv (p 2 ; A uv ) = 0. (b) The UV-cut SD kernel K vv (p 2 ; Ag$) 
for the critical case on DCSB. The dotted curve denotes the 
original SD kernel function K{p 2 ). 
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FIG. 5: (a) The infrared quark mass M m (0; Air)/M in the 
IR-cut SD equation plotted against the artificial IR-cutoff pa- 
rameter Air. For Air <0.4GeV, there is no significant effect 
observed for DCSB. For Air > A m !t ~ 0.53GeV, no DCSB 
is observed as Mi R (p 2 ;Ai R ) = 0. (b) The IR-cut SD kernel 
K m (p 2 ; Ajr') for the critical case on DCSB. The dotted curve 
denotes the original SD kernel function K(p 2 ). 
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FIG. 6: (a) The infrared quark mass Mim(0;c)/Mo in the 
IM-suppressed SD equation plotted against the artificial IM- 
suppression parameter c. No DCSB is found as Mim(p 2 ; c) = 
for c < c crit ^ 0.58. (b) The IM-suppressed SD kernel 
Kim(p 2 ; Ccrit) for the critical case on DCSB. The dotted curve 
denotes the original SD kernel function K(p 2 ). 




FIG. 7: The infrared quark mass M(0; A) in the scaled SD 
equation against A. For A > 0.8, the quark mass function 
M(p 2 ; A) is an increasing function on A for fixed value of p 2 . 
For A < A crit ^ 0.8, no DCSB occurs as M(p 2 ; A) = 0. 
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FIG. 8: The thermal infrared quark mass M T (p 2 = 0,lo'q) 
obtained from the thermal SD equation plotted against the 
temperature T. Chiral symmetry restoration is found at a 
critical temperature T c ~ lOOMeV. 
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FIG. 9: The thermal quark mass function Mt(p 2 ,u^) at 
various Matsubara frequencies n = 0,1,2,3 obtained from 
the thermal SD equation for (a) a low temperature case of 
T = 0.6T C = 60MeV and (b) a high temperature case of 
T = 0.9T C = 90MeV. 
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FIG. 10: The comparison of the thermal quark mass function 
M t (p 2 ,loI) with M T (p 2 ) = Mt(p 2 + u„) for n = 1,2,3 at 
(a) a low temperature T = 0.6T C = 60MeV and (b) a high 
temperature T = 0.9T C = 90MeV. For each n, the dashed 
curve denotes Mt(p 2 ) defined from M(p 2 ,u> 2 ). 
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FIG. 11: The thermal quark condensate (qq(T))A/{qq(0))A 
plotted against the temperature T for A = 4GeV. 



